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Aircraft noise mitigation is an ongoing challenge for the aeronautics research community. 
In response to this challenge, low-noise aircraft concepts have been developed that exhibit 
situations where the jet exhaust interacts with an airframe surface. Jet flows interacting 
with nearby surfaces manifest a complex behavior in which acoustic and aerodynamic char- 
acteristics are altered. In this paper, the variation of the aerodynamics, acoustic source, 
and far-field acoustic intensity are examined as a large flat plate is positioned relative to 
the nozzle exit. Steady Reynolds- Averaged Navier-Stokes solutions are examined to study 
the aerodynamic changes in the field-variables and turbulence statistics. The mixing noise 
model of Tam and Auriault is used to predict the noise produced by the jet. To validate 
both the aerodynamic and the noise prediction models, results are compared with Par- 
ticle Image Velocimetry (PIV) and free-field acoustic data respectively. The variation of 
the aerodynamic quantities and noise source are examined by comparing predictions from 
various jet and flat plate configurations with an isolated jet. To quantify the propulsion air- 
frame aeroacoustic installation effects on the aerodynamic noise source, a non-dimensional 
number is formed that contains the flow-conditions and airframe installation parameters. 


Nomenclature 


Abbreviations 

CFD 

Computational Fluid Dynamics 

FUN3D 

Fully Unstructured Navier-Stokes Three-Dimensional Solver 

HWB 

Hybrid Wing Body 

NPR 

Nozzle Pressure Ratio 

OASPL 

Overall Sound Pressure Level 

PAA 

Propulsion Airframe Aeroacoustics 

PIV 

Particle Image Velocimetry 

PSD 

Power Spectral Density 

RANS 

Reynolds- Averaged Navier-Stokes 

SHJAR 

Small Hot Jet Acoustic Rig 

SMC 

Small Metal Chevron 

SPL 

Sound Pressure Level 

SST 

Shear Stress Transport 

TICE 

Turbulent Kinetic Energy 

TTR 

Total Temperature Ratio 

Symbols 

A 

Empirical constant 

c 

Speed of sound 

Coo 

Free-stream speed of sound 
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Empirical constant associated with turbulent time scale 
Empirical constant associated with turbulent length scale 
Nozzle exit diameter 
Fully expanded jet diameter 

Turbulent kinetic energy associated with fine-scale mixing noise 

Turbulent length scale 

Acoustic Mach number 

Design Mach number 

Fully expanded Mach number 

Free-stream Mach number 

Static pressure 

Acoustic pressure 

Fully expanded static pressure 

Stagnation pressure 

Statistical source term 

Strouhal number 

Spectral density of acoustic pressure 
Stagnation temperature 
Fully expanded static temperature 
Mean streamwise velocity component 
Fully expanded jet velocity 
Vector observer position 
Estimate of jet impingement location 
Axial distance to the plate trailing edge 
Source vector 

Radial distance from the plate to the jet centerline 

Spreading angle of the jet 

Dissipation of turbulent kinetic energy 

Non-dimensional number 

Ratio of specific heats 

Turbulent time scale 

Radial frequency 

Observer angle from the upstream axis 


Introduction 

Aircraft noise mitigation is an ongoing challenge for the aeronautics research community. Noise associated 
issues such as acoustic fatigue, annoyance to the public, and hearing loss have led to increased efforts toward 
the reduction of jet noise. Increased understanding of jet noise produced in various flow conditions allows 
for improved acoustic prediction models, that are essential for the development of new low-noise aircraft 
designs. In support of these aircraft designs, shielding techniques have been studied that involve airframe 
structures located near jet engines to minimize noise perceived by the observer on the ground (see for example 
Papamoschou and Mayoral 1 or Huang and Papamoschou 2 * ). Jet flows interacting with these nearby surfaces 
result in complex scattered acoustic fields that are difficult to model. 

Supersonic aircraft configurations, such as those developed by Welge et al A 4 for the 2020 to 2035 
timeframe or by Morgenstern et al. 5 for the 2035 timeframe and beyond, demonstrate complicated nozzle 
and engine placement designed to reduce sonic boom and jet noise. Advanced low bypass ratio engines, 
as shown by Sokhey and Kube-McDowell, 6 can be integrated into the airframe while potentially reducing 
jet noise. These concepts could also exhibit situations where the jet exhaust interacts directly with an 
airframe surface resulting in additional noise sources and a complicated scattered field. New jet noise 
prediction methodologies need to be developed to address these complicated configurations and advanced 
nozzle designs. 

Another advanced concept is the Hybrid Wing Body (HWB) aircraft that shields jet noise by placing the 
engines above the airframe. This effect is described by Thomas et al. ' and Czech et al. 8 Several techniques 
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have been developed to predict the scattered noise from the engine. Some of these techniques include the 
Boundary Element Method, 2 the Equivalent Source Method used in NASA’s Fast Scattering Code, 9 and the 
ray tracing method. 10 However, only a limited number of source models have been developed for use with 
scattering methodologies. 

Jet surface interactions produce a change in the noise sources associated with the aerodynamics of the 
jet flow. There have been several studies that examined the aerodynamic effects of nearby and impinging 
surfaces on jet flows. Donaldson and Snedeker 11 and Lamont and Hunt 12 showed that jets impinging on 
oblique and perpendicular surfaces significantly affect the turbulence and shock cell structure of supersonic 
jets. Studies involving jets not directly impinging on surfaces include those by Sawyer 13 and Al-Qutub and 
Budair 14 among others. Capturing the effect on the acoustic source due to the changes in the aerodynamic 
characteristics of the jet flow-field is essential for accurate prediction of jet noise. The effect of nearby surfaces 
on the aerodynamic flow-field is dependent on multiple parameters such as NPR, TTR, nozzle geometry, the 
boundary condition on the surface, and the position of the surface relative to the jet. 

In this paper, the variation of the aerodynamics, acoustic source, and far-field acoustic intensity are 
examined as a large flat plate is moved relative to the nozzle exit. Steady Reynolds-Averaged Navier-Stokes 
(RANS) solutions are found using realistic nozzle and flat plate geometries to study the aerodynamic changes 
in the field-variables and turbulence statistics. The semi-empirical fine-scale mixing noise model of Tam and 
Auriault 15 is used to predict the noise produced by the jet. This model is also used to examine the integrand 
of the acoustic source that uses quantities obtained from the steady RANS solution as its argument. To 
validate both the aerodynamic and the noise models, results are compared with Particle Image Velocimetry 
(PIV) and free-field acoustic data, respectively. The variation of the field-variables, turbulence statistics, 
and noise is examined by comparing predictions from various jet and plate configurations with an isolated 
jet. Multiple jet operating conditions are also examined. A non-dimensional number is developed to quantify 
the effect of the surface on the aerodynamic noise source based on the jet operating condition and surface 
location. 

In the following sections of this paper, a non-dimensional number is introduced to quantify the effect of the 
airframe on the aerodynamic source, and the associated parameters describing the jet condition and surface 
location are defined. Next, the steady RANS model is presented and compared with the PIV data of Bridges 
and Wernet. 16 Solutions of the steady RANS model are shown for multiple jet and plate configurations. 
Following the steady RANS assessment, a statistical noise prediction approach is described and compared 
with measurements of Bridges and Brown. 1 ' Then, noise predictions of the isolated jet are compared to 
the jet with a flat plate at various positions. Finally, a conclusion of this study and related future work are 
discussed. 


Analysis 

The method by which a jet engine is integrated with the airframe of a flight vehicle can have significant 
effect on the aerodynamic source of sound. This is a highly complicated problem that has received significant 
attention. 7,8, 18 It is characterized by a large number of parameters that are highly interdependent. A non- 
dimensional number is formed, which has arguments involving the flow-conditions and jet position, that can 
be used as a basic guide to ascertain whether the aerodynamic source is affected by the airframe relative to 
the equivalent isolated jet aerodynamic source. 

The coordinate system and geometry are illustrated in Fig. 1. Coordinates x, y, and z are normalized by 
the nozzle diameter. This normalized coordinate system is used to illustrate results in the following sections. 
The origin of the coordinate system used for the analysis is located at the nozzle exit. The positive x axis 
points in the jet principal direction, the y axis is normal to the plate, and the z axis is parallel to the plate. 
The plate is located at multiple cross-stream distances from the jet centerline axis to the plate surface, y p , 
and multiple streamwise distances from the nozzle exit to the trailing edge, x p . 

Parameters are identified that have a direct impact on the aerodynamic source of the jet based on the 
position of the airframe surface. These parameters are illustrated in Fig. 1. The first parameter is the fully 
expanded diameter of the jet, 
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where D is the nozzle exit diameter, Dj represents the necessary equivalent nozzle exit diameter for a shock 
free flow, M r j and Mj represent the design Mach number and the fully expanded Mach number, respectively, 
and 7 is the ratio of specific heats. 

Additionally, x p is the distance from the nozzle exit to the trailing edge of the airframe 3, parallel to the 
jet centerline. y v is the characteristic length from the nozzle centerline to the nearest airframe surface, xj 
is the shadow distance on the airframe surface from the nearest airframe surface point at which the jet first 
interacts with the airframe. This value is dependent on the problem geometry, the engine cycle condition, and 
Xi . The simplest estimate for xi is the jet impingement location, that is a function of the initial spreading 
rate of the jet. Here xi is approximated as, 


_ y P - -P/2 

tan [(>,,] 


(2) 


where S v is the spreading angle of the jet. Values of S v are not readily available without numerical calculations 
or measurement. Here, we adopt an empirical model developed by Lau 19 for S v , 


S v = 0.177(1 - 0.294 M 2 ) ^1 + ^(M 2 - 1 ){Tj/T 0 - 1.4) 2 ^ , (3) 

which is valid for a wide range of single-stream jet Mach numbers and temperature ratios. The flows of 
this investigation fall within the range of validity of Eqn. 3. Equation 3 is dependent on a wide range of 
parameters that do not explicitly appear. Using these parameters a non-dimensional number is proposed, 



Alternatively, Eqn. 4 can be written as, 


(4) 


_ Dj Xp tan[(5^] 

~y P {y P -D/2Y 

where y p > D/2. Physically, Eqn. 5 is the ratio of the product of the jet and airframe length scales divided 
by the cross-stream length scale and interaction distance. Small values of T can suggest that airframe effects 
on the jet aerodynamic noise sources are negligible. Likewise, large values of T can suggest that the effects 
of the airframe on the aerodynamic noise sources of the jet are very large. In the following sections, Eqn. 5 
is evaluated using different jet conditions and airframe surface positions, and compared with numerical 
predictions of the variation of the noise from the aerodynamic source. 


Results 


Steady RANS Assessment 

The steady RANS equations are solved with a Computational Fluid Dynamics (CFD) approach using the 
NASA Fully Unstructured Navier-Stokes three-dimensional 20 (FUN3D) solver. The Menter 21 Shear Stress 
Transport (SST) turbulence model is used to close the RANS equations. The model utilizes the strengths 
of the Jones and Launder k-e model 22 and the Wilcox k-u> model . 23,24 It utilizes the k-iv model in the inner 
region of the boundary layer and switches to the k-e model in the outer region and in free shear layers. The 
model also takes into account the transport of the principal turbulent shear stress within adverse pressure 
gradients. In this work, the noise source is dependent upon quantities obtained from the steady RANS 
solution. 

The plate positions relative to the nozzle exit assessed in this paper are shown in Table 1. The first 
column contains normalized streamwise distances from the nozzle exit to the edge of the flat plate and the 
second column shows normalized radial distances to the flat plate. Table 2 provides information regarding 
the jet operating conditions analyzed. The operating conditions include three subsonic jets, one of which is 
heated, and two cold supersonic jets, one of which is over-expanded and the other ideally-expanded. 

The field variables and turbulent statistics are discretized on a mixed element (structured-unstructured) 
computational grid. An unstructured grid is used to resolve the jet flow when the plate is present while a 

a We elect to use airframe and flat plate interchangeably in terms of the analysis. 
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structured grid is retained to resolve the jet plume. All calculations are fully three-dimensional and utilize 
symmetry. The nozzle geometries and flow conditions coincide with the PIV experiment of Bridges and 
Wernet. 16 Figure 2 shows an outline of the three-dimensional computational domain with the plate located 
two nozzle diameters, D, offset from the jet axis and parallel to the jet centerline. Values of p Q and T 0 are 
specified by boundary conditions consisting of nozzle pressure and temperature ratios at the nozzle inlet 
and static pressure at the domain exit. A symmetric boundary condition is enforced on the x-y plane at 
z = 0. All other boundaries are defined with a free-stream condition of =0.01 and ambient pressure. 
The SMC000 nozzle profile and a portion of the computational domain are shown in Fig. 3. This portion of 
the domain represents a slice in the x-y plane at z = 0 normalized by D. The full computational domain 
extends 75 D downstream, 501? cross-stream along the z axis, and 100Z? cross-stream along the y axis to 
form a rectangular prism. The number of grid points in the domain with a plate present is 2,195,937 and 
the number of elements is 2,827,906. The number of grid points in the domain without a plate present is 
1,180,575 and the number of elements is 1,472,843. These represent typical values for the cases examined. 

Validation of the isolated jet steady RANS solution is performed by comparing the streamwise velocity 
component and Turbulent Kinetic Energy (TKE) with the PIV data from the isolated jet experiment of 
Bridges and Wernet. 16 The PIV dataset of Bridges and Wernet 16 was developed in part to support the 
validation of jet noise prediction methodologies and steady RANS solutions for a wide range of jet velocities 
and temperatures. Mean velocities, TKE, and Reynolds stresses are important parameters for statistical 
acoustic analogies. The first two are used for validation and subsequent analysis in this study. 

The first comparison of steady RANS solutions with measurement is performed with a subsonic cold jet 
operating at Mj = 0.513 from the convergent SMC000 nozzle. Contour plots of Fig. 4 show a qualitative 
comparison of the streamwise velocity component and TKE. The velocity and TKE are normalized by the 
fully expanded velocity and its square, respectively. The steady RANS solutions are shown on the top half 
planes and the corresponding PIV data are shown on the bottom half planes. The streamwise velocity shown 
in Fig. 4(a) has noticeable variation from the PIV. The predicted thickness of the potential core past 6 D 
is below measurement and the fall-off past 8 D is higher. The predicted peak TKE occurs near 5.5 D and 
the PIV peak TKE occurs further downstream at 6.251?. Furthermore, the predicted magnitude of peak 
TKE is lower than measurement. However, the global qualitative agreement is favorable compared to other 
solutions produced by similar CFD codes (for example see Georgiadis et al. 25 ). 

Qualitative comparisons are conducted for the isolated jet by examining aerodynamic data along the 
centerline and axial locations at x/ D = 1, 4, and 16. Figure 5 shows the centerline variation of the normalized 
streamwise velocity component and TKE of both the steady RANS solution and PIV data. Figure 6 shows 
radial profile comparisons of the same quantities. Figure 5(a) shows that the jet potential core length is 
larger than measurement by 2D. The predicted streamwise velocity component follows the same ss \/x 
decay as the measured data along the centerline, shown in Fig. 5(a), and also matches the radial decay at 
x/D = 1 and x/D = 4, shown in Fig. 6(a). The predicted TKE along the centerline matches the experiment 
in terms of peak magnitude but the peak occurs one diameter further upstream than experiment. The peak 
magnitudes of TKE at x/D = 1 and x/D = 4 are larger than measurement by 0.0015 and 0.001 TKE /U 2 
respectively and are located 0.251? closer to the centerline than experiment. The solution under-predicts 
both the streamwise velocity component and TKE relative to measurement far downstream from the nozzle 
exit. 

Qualitative comparisons are also conducted for the supersonic over-expanded jet operating at Mj = 1.29 
from the M r j = 1.5 SMC016 nozzle. Again, the aerodynamic data is examined along the centerline and 
radial locations at x/D = 1,4, and 16. Figure 7 shows the centerline variation of the normalized streamwise 
velocity component and TKE of both the steady RANS solution and PIV data. Figure 8 shows radial profile 
comparisons of the same quantities. As seen in Fig. 7, the predicted potential core length is 2D larger relative 
to measurement. Within the potential core, the predicted shock cell structure locations compare favorably 
with the experiment. The steady RANS solution over-predicts the rate of decay along the centerline shown 
in Fig. 7(a) and matches the radial decay in the potential core region at x/D = 1 and 4 as shown in Fig. 8(a). 
The steady RANS solution under-predicts the radial decay through the shear layer greater than 0.5Z? from 
the centerline. In the fully developed region of the flow at x/D = 16, the steady RANS solution only slightly 
over-predicts the radial velocity profile. This is a favorable prediction relative to the TKE profile at x/D = 16 
shown in Fig. 8(b) and the subsonic jet predictions shown at x/D = 16 in Fig. 6. On the centerline the 
peak TKE/Uj 2 is 0.03 below measurement and occurs 2D further downstream than measurement as shown 
in Fig. 7(b). The TKE predictions match the experiment in peak magnitude at x/D = 4 but over-predicts 
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TKE/Uj 2 at x/D = 1 and x/D = 16 by 0.015 and 0.005, respectively. In concurrence with the subsonic cold 
jet comparisons of Fig. 6(b), the predicted TKE is lower than measurement further downstream from the 
nozzle. These trends are representative of all jet conditions shown in Table 2. The predictions are favorable 
relative to experiment, for all jets examined, in the potential core region of the flow-field where inviscicl terms 
dominate the equations of motion. 

The steady RANS solutions of a Mj = 0.513 and TTR = 1.00 jet from the convergent SMC000 nozzle 
in the presence of a flat plate are now compared with the isolated jet operating at the same conditions. 
Figure 9 shows radial profiles of TKE at x/D = 10 for multiple plate positions. The plate is located at 
y p /D = — 1, — 2, — 4, — 6, — 8, and —10 perpendicular to the jet centerline and extends to x p /D = 10 and 
20 downstream from the nozzle exit. Figure 9 shows that as the plate is moved closer to the jet the TKE 
distribution is increasingly deformed. The peak magnitude closest to the plate decreases and the peak 
magnitude furthest from the plate increases. This trend is amplified when increasing the x p /D location from 
10 to 20 for most of the cases examined. Furthermore, the jet plume is being deformed and is drawn toward 
the plate due to a coanda like effect. For example, in Fig. 9(a) the x v /D = 10 and isolated jet case both 
have a TKE local minimum at y/D = 0. However, as the plate is moved closer as in Fig. 9(f) the local 
minimum of TKE for x p /D = 10 is now located at y/D = —0.16. The effect of extending the plate from 
x p /D = 10 to Xp/D = 20 amplifies the effect of the deformation of the jet plume as can be seen in Fig. 9(a). 
The peak TKE close to the plate is lower by 0.0005 TKE/Uj 2 , and the peak TKE furthest from the plate 
is higher by 0.001 TKE/Uj 2 . Similar changes of the TKE distribution are observed with all jet conditions 
examined. These preliminary numerical results show the changes in the aerodynamic characteristics of the 
jet plume induced by nearby surfaces. These changes effect the noise source amplitude and position as the 
plate location is altered. 


Aeroacoustic Assessment 


The results of the aerodynamic assessment show that nearby solid surfaces change the aerodynamic charac- 
teristics of the jet flow-field, even if the jet centerline is multiple diameters away from a solid surface. To 
assess the changes of the acoustic source quantitatively, a statistical noise prediction approach is selected. 
Here, we choose the semi-empirical fine-scale mixing noise model of Tam and Auriault. 15 The expression for 
the spectral density of the acoustic pressure in the far-held, S, is given by, 


S< a! , u )= 4 »(j^) 


3/2 /.oo r- oo /.oo ~2|3 \p a (x\ J/,w)| 2 exp 
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where A is a constant coefficient, c is the speed of sound, u is the mean streamwise velocity component, x 
is the observer position, y is the source position, and u> is the radial frequency. q s = (4/9 )A 2 c 2 ~p 2 k 2 is a 
statistical source term where k s is the TKE associated with turbulence that produces fine-scale mixing noise. 
k s is set equal to the TKE computed from the steady RANS solution. 

If predictions are restricted to the sideline direction and it is assumed that sound refraction by the jet 
shear layer has negligible effect on S, then a simplified form of p a { x \ y ■> w) can be constructed. As shown by 
Morris and Farassat 26 the adjoint acoustic pressure at ip = 7r/2 (the jet sideline) is, 


\Pa{X2',y,u) | 2 


UJ 


64-7T 4 




( 7 ) 


This approach does not take into account propagation effects associated with sound waves interacting with 
surfaces such as the flat plate. This is advantageous because the effects of surfaces on the jet noise source 
are isolated. 

The scales of turbulence in Eqn. 6 are related to the steady R ANS solution of the jet by simple dimensional 
models. The length scale is l s = cikV 2 /e and the time scale is r s = c T fc s /e, where e is the dissipation of 
TKE. The quantity, e, is computed from the steady RANS solution as e = 0.09 k s ui. By relating the turbulent 
scales to the steady RANS solution, the predicted noise is dependent on the jet mean flow. The mean flow 
is dependent on the boundary conditions, the nozzle geometry, and the position of the plate relative to the 
jet. The acoustic source strength and spatial distribution as a function of frequency is directly connected to 
the steady RANS solutions. 
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The coefficients A, ci , and c T of Eqn. 6 are specified by comparing prediction with measurement, following 
the methodology of Tam and Auriault. 15 The coefficients are calibrated with the SMC000 nozzle operating 
at Mj = 1.00 and TTR = 1.00. The coefficients are constant irrespective of observer position, jet operating 
condition, or airframe geometry after calibration. The values are A = 3548.0, C; = 0.018, and c T = 0.015, that 
vary from those calculated by Tam and Auriault 15 because the steady RANS solver and turbulence model 
differ. Figure 10 shows the prediction compared with the experiment of Bridges and Brown. 1 ' The Sound 
Pressure Level (SPL) per unit St is represented on the y-axis, where St = ujDj/uj is the Strouhal number. 
The prediction captures the peak intensity and the nearly correct decay of intensity at high frequency. There 
is slight discrepancy at lower frequencies where the decay of intensity is lower than measurement. 

The effect that the flat plate has on acoustic intensity originating from the jet aerodynamic source is 
investigated by using the steady RANS solutions and Eqn. 6. The isolated jet predictions are subtracted from 
the installed jet predictions on a Power Spectral Density (PSD) basis for comparison. The first comparison 
consists of the cold subsonic jet at Mj = 0.513 and is presented in Fig. 11. In Fig. 11(a), the plate is 
located at y p /D = —1, —2, —4, —6, —8, and —10 from the jet centerline and extends x p /D = 10 downstream 
of the nozzle exit. In Fig. 11(b), the plate is located at y p /D = —1 from the jet centerline and extends 
x p /D = 4, 10, and 20 downstream of the nozzle exit. When the plate is extended to x p /D = 10 downstream 
of the nozzle exit, the results from Fig. 11(a) show only a small effect on the noise spectrum as the plate 
approaches the centerline until the plate is at y p /D = —2. The noise deviation from the isolated jet reaches a 
maximum of —1.5 dB at the lowest frequency and 0.75 dB at the highest frequency when the plate is located 
at y p /D = —2. For all plate locations further than 2D away from the centerline, the noise deviations are 
within ±0.75dB from the isolated jet case. When the plate is located at y p / D = —1, the noise intensity is 
—6.9 dB relative to the isolated jet case at the lowest frequency. The effect of varying the plate length relative 
to the nozzle exit is displayed in Fig. 11(b). At the lowest frequency, the noise intensity deviation from the 
isolated jet reaches a maximum of —6.5 dB, —6.9 dB, and —7.9 dB as the plate is extended to x p /D = 4, 10, 
and 20 downstream from the nozzle exit, respectively. At the highest frequency, the maximum deviation for 
each plate extension does not exceed 0.75 dB. 

Predictions of the over-expanded supersonic Mj = 1.29 jet from the convergent-divergent SMC016 nozzle 
are compared at various plate positions relative to the isolated jet. Figure 12(a) displays the effects of the 
surface as it approaches the jet centerline, and Fig. 12(b) shows the effects of the surface as the plate length is 
extended relative to the nozzle exit. In Fig. 12(a), the noise intensity deviation from the isolated jet reaches a 
maximum of —1.0 dB at the lowest frequency and OAdB at the highest frequency when the plate is located at 
y p /D = —2. For all plate locations further than 2D away from the centerline, the noise deviations are within 
—0.6 dB from the isolated jet case at the lowest frequency and — 0.02dB at the highest frequency. When 
the plate is located at y p /D = —1, the noise reaches a maximum —OAdB difference from the isolated jet 
case. From Fig. 12(b), the noise deviation from the isolated jet reaches a maximum of —5.7 dB, —OAdB, and 
—7.7 dB at the lowest frequency as the plate extends x p /D = 4, 10, and 20 downstream from the nozzle exit, 
respectively. The maximum dB difference from the isolated jet for each plate position is slightly decreased 
in magnitude at lower frequencies when compared to the subsonic jet. There is no increase in noise intensity 
at high frequencies in contrast to the subsonic case. 

Fig. 13 shows comparisons of an on-design Mj = 1.5 cold jet. From Fig. 13(a), the noise intensity 
deviation from the isolated jet reaches a maximum of —1.0 dB at the lowest frequency and OAdB at the 
highest frequency when the plate is located at y p /D = —2 and x p /D = 10. For all plate locations further 
than 2D from the centerline, the noise deviations are within —0.6 dB from the isolated jet case at the lowest 
frequency and — 0.02dB at the highest frequency. When the plate is located at y p /D = —1, the noise reaches 
a maximum — OAdB difference from the isolated jet case at the lowest frequency. In Fig. 13(b), the noise 
deviation from the isolated jet reaches a maximum of — OAdB, — OAdB , and — 8AdB as the plate extends 
to x p /D = 4, 10, and 20 from the nozzle exit, respectively. In regards to the magnitude of the maximum 
difference and general trend of the spectra, the predictions of the jet and plate cases relative to the isolated 
jet case are in close agreement with the over-expanded jet. In Fig. 13(b), there is a OAdB increase in 
magnitude of noise intensity deviation with the plate located at y p /D = —1 and extending x p /D = 4 and 
20 downstream when compared to the over-expanded jet in Fig. 12(b) at lower frequencies. At the higher 
frequencies there is no notable difference from the over-expanded jet. 

The comparisons described above are representative of the other jet conditions analyzed in this study. As 
the jet approaches the plate there is a consistently larger difference in the noise spectrum from the isolated 
jet at lower frequencies. This is due to the surface having a larger effect on the flow further downstream 
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from the nozzle exit. It has been shown by Brooks et al. 2 ‘ and Podboy 28 that peak noise sources are located 
near the nozzle exit at higher frequencies and lower frequency peak noise sources are located multiple nozzle 
diameters downstream. The presence of the plate has a larger effect on the jet flow multiple diameters from 
the nozzle exit and therefore has a larger effect on the aerodynamic source at lower frequencies. 

Next, the effect that the plate has on the acoustic source localization is examined by comparing noise 
source maps of the jet flow for various plate positions using the steady RANS solutions and the integrand 
of Eqn. 6 near the peak frequency. Contours of SPL per unit St for the cold subsonic Mj = 0.513 jet 
at a frequency of 1 kHz (St ss 0.3) are shown in Fig. 14. In this comparison the plate length extends to 
Xp/D = 10 and the plate is located at y p /D = —10, —6, —4, —2, and —1 from the jet centerline. When the 
plate is located at y p /D = —10 in Fig. 14(a), the noise source distribution is unaffected by the presence of 
the plate. The two peak noise sources are located 9.5 D downstream from the nozzle exit and are symmetric 
in magnitude about the jet centerline. There is no significant deformation of the source distribution due 
to the presence of the plate further than 2D from the jet centerline. As the plate approaches the jet at 
y p /D = —2 shown in Fig. 14(d), the peak magnitude of the source closest to the plate is decreased by ~ 1 cLB 
and the location is unchanged. When the plate is located at y p /D = —1, shown in Fig. 14(e), the peak noise 
source furthest from the plate is decreased in magnitude by ss ldB and is shifted 0.51? upstream relative 
to the y p /D = —10 case shown in Fig. 14(a). The peak noise source distribution closest to the plate is 
significantly deformed. The peak magnitude is decreased by « 6 dB and is shifted ID upstream relative to 
the y p /D = — 10 case. An additional strong noise source is also formed close to the plate further downstream 
at x/D = 14. This is a distinguishable difference from the other comparisons and could be explained by an 
observed increase in TKE past the trailing edge of the plate. Figure 14 shows that the magnitude of the 
peak acoustic source closest to the plate decreases as the plate surface approaches the jet centerline, and the 
magnitude of the peak acoustic source furthest from the plate is not significantly altered until the plate is 
located at y p /D = —1 from the jet centerline. The location of the peak noise sources are unaltered until the 
plate is located at y p / D = — 1. These comparisons represent trends that are consistent for all jet conditions 
analyzed. 

Overall, the results show that the plate has a larger effect on the acoustic source as it approaches the 
jet centerline and as the plate length extends further downstream. As the effective jet impingement area of 
the plate is increased, the acoustic intensity radiating from the jet decreases. This result does not account 
for additional sources produced by the jet interacting with the surface. To quantify the effect of a nearby 
surface on the jet noise source, the non-dimensional number T is used, as described in Eqn. 5. The overall 
sound pressure level (OASPL) is predicted over a frequency range of 20 Hz to 100 kHz for all plate locations 
and jet conditions described in Tables 1 and 2 using the steady RANS solutions and Eqn. 6 at ip = 7r/2. The 
installed jet predictions of OASPL are subtracted from the isolated jet predictions on a PSD basis. Figure 15 
shows the absolute value of AOASPL as a function of T. The physics of the jet flow and airframe interaction 
dictate the value of T, and the value of T expresses how large an effect the jet airframe interaction has on 
the jet noise source. The data collapse shows a critical value of T « 1. As the cross-stream length scales 
VpiUp ~ -D/2) decrease to the equivalent of the product of the jet and airframe length scales (DjX p tan[d^]), V 
increases from 0 to 1. In the region, 0 < Y < 1, the cross-stream length scales are dominant and the physical 
quantities of the aerodynamic flow are not sufficiently altered to have a significant effect on the jet noise 
source. In Fig. 15, small values of T approaching ss 1 have AOASPL less than 0.5 dB and are considered 
negligible. As the jet spreading angle, fully expanded jet diameter, or surface length increases the product 
of the jet and airframe length scales surpasses the cross-stream length scales and T increases beyond unity. 
As a result of the dominance of the jet and airframe length scales in this region, 1 < T < oo, the physical 
quantities of the aerodynamic flow are sufficiently altered and have a large effect on the jet noise source. As 
shown in Fig. 15, values of T > 1 result in a range of AOASPL from 2.0 dB to 2.75 dB. It is observed that 
AOASPL increases as T increases. 

Calculation of T can be a useful tool when performing propulsion airframe aeroacoustics (PAA) analysis. 
When analyzing the noise propagation for a jet airframe interaction corresponding to T « 1, the airframe 
is expected to have a negligible effect on the jet noise source. For instance, a supersonic on-design Mj = 1.5 
cold jet with an airframe surface located at x p /D = 20 and y p /D = 4 corresponds to T = 0.099 and AOASPL 
= 0.083dB in Fig. 15. Since T << 1 and there is a negligible difference in the free-fielcl jet noise source, the 
isolated jet aerodynamic source model along with a tailored Green’s function can used for a prediction that 
includes propagation about the airframe. A M 3 = 0.985 jet with an airframe surface located at x p /D = 4 
and y p /D = 1 corresponds to T = 5.063 and AOASPL = 2A7dB in Fig. 15. The large value of T results 
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in a large difference from the free-stream noise prediction; therefore a separate aerodynamic source model 
is required for further analysis. An isolated jet aerodynamic source model can be used for a PAA analysis 
with small values of T. For large values of T, some noise reduction can be attributed to the change in the jet 
noise source and not to shielding effects. Recent studies have investigated the effects of liners on propulsion 
airframe surfaces. 18 The effectiveness of airframe liners is associated with the location of the airframe relative 
to the jet flow. The implications of T could be useful in future PAA liner studies. 

T can include additional terms involving multiple jet streams or account for the boundary layer thickness 
on the airframe. An equivalent parameter for the fully expanded jet diameter and jet spreading angle can 
be developed for a dual-stream jet. However, a new model would be required for the estimation of the jet 
spreading angle. Inclusion of boundary layer effects will increase the effect of the cross-stream length scales. 
The airframe boundary layer thickness can be subtracted from the cross-stream length scales. A model to 
estimate the boundary layer thickness would be required. The general effect of the airframe on the jet noise 
source as a function of T is expected to be consistent for these cases but may have a different critical value. 

Conclusion 

This study is motivated by low-noise aircraft designs that exhibit situations where the jet exhaust interacts 
with an airframe surface and to develop an understanding of the effect of nearby surfaces on the aerodynamic 
noise sources of jet flows. A simplified analysis is performed to examine the variation of the aerodynamics, 
acoustic source, and far-field acoustic intensity as a large flat plate is moved relative to the nozzle exit. 
Steady RANS solutions are found using laboratory nozzle and flat plate geometries to study the aerodynamic 
changes in the field- variables and turbulence statistics. The semi-empirical fine-scale mixing noise model of 
Tam and Auriault is used to predict the noise produced by the jet and to examine the integrand of the 
acoustic source. The aerodynamic and noise prediction models are validated by comparing results with PIV 
and free-fielcl acoustic data, respectively. The variation of the field-variables, turbulence statistics, and noise 
are examined by comparing predictions from various plate configurations with an isolated jet. The analysis 
includes multiple jet operating conditions. 

Results of the aerodynamic assessment show that nearby solid surfaces change the aerodynamic char- 
acteristics of the jet flow-field, even if the jet centerline is multiple diameters away from a solid surface. 
The effect of the surface is amplified as it approaches the jet centerline and as the plate length extends 
further downstream. The induced change in the aerodynamic flow-field is shown to have a direct effect on 
the source of jet noise. The installed noise predictions relative to the isolated jet show that the acoustic 
intensity originating from the jet aerodynamic source is decreased as the effective impingement surface area 
of the plate is increased. 

A non-dimensional number is developed to quantify the effect of the surface on the jet aerodynamic source 
based on the jet condition and surface location parameters. The reference number T physically describes the 
ratio of the product of the jet and airframe length scales over the cross-stream length scales. The evaluation 
of T over the range of surface locations and jet conditions yields a critical value of T « 1. Small values of 
T suggest that airframe effects on the jet flow are negligible. Large values of T imply that the effects of the 
surface on the jet aerodynamic source are very large. This simple calculation of T can be used as a basic 
guide to determine if the aerodynamic source is affected by the airframe relative to the equivalent isolated 
jet aerodynamic source. Additional terms can be included with T to capture additional physical effects. 
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Tables 


Table 1. Plate Locations 


Axial Location of Trailing Edge ( x p /D ) 

Radial Location Relative to Jet Centerline ( y p /D ) 

4 

1,2,4,6,8,10 

10 

1,2,4,6,8,10 

20 

1,2,4,6,8,10 


Table 2. Jet Operating Conditions 


Nozzle 

Setpoint 

NPR 

TTR 

Ma 

Mi 

SMC000 

3 

1.197 

1.000 

0.5 

0.513 

SMC000 

7 

1.861 

1.000 

0.9 

0.985 

SMC000 

27 

1.357 

1.814 

0.9 

0.678 

SMC016 

11606 

2.748 

0.761 

1.128 

1.29 

SMC016 

11610 

3.670 

0.706 

1.31 

1.5 


Figures 


y 

4 



Figure 1. The coordinate system and distance parameters used throughout this study. 
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Figure 2. An outline of the computational domain used for the CFD calculations. The edges of the compu- 
tational domain are shown as black lines. The flow of the jet is in the positive x direction and the coordinate 
system is centered at the nozzle exit. The domain extends 75 nozzle diameters downstream, 100 diameters in 
the y direction, and 50 diameters in the z direction. The types of boundaries are labeled. 



-5 - 4 - 3 - 2-10 1 2 3 


x/D 


Figure 3. A portion of the computational domain on the x-y plane at z = 0 showing the structured/unstructured 
grid of the convergent nozzle and flat plate located 2D away from the centerline. The nozzle geometry matches 
the NASA Glenn Research Center SMC000 nozzle with an exit diameter of D = 0.0508 m. The coordinates are 
normalized by the nozzle exit diameter. 
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(a) Streamwise Component of Velocity 
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(b) Turbulent Kinetic Energy 


Figure 4. Contours of (a) the streamwise velocity component and (b) TKE of the steady RANS solution and 
PIV data from the experiment of Bridges and Wernet. 16 The steady RANS solutions are shown on the top 
half planes and PIV data are shown on the bottom half planes. The streamwise velocity component and TKE 
are normalized by the fully expanded velocity and the square of the fully expanded velocity respectively. The 
coordinates are normalized by the nozzle exit diameter D = 0.0508 m. The jet operates at Mj = 0.513 and 
TTR = 1.00 from the convergent SMC000 nozzle. 
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Fun3D 




(a) Velocity (b) TKE 


Figure 5. The (a) steady RANS streamwise velocity component and (b) TKE along the jet centerline (y/D = 0 
&; z/D = 0) compared with PIV data from the experiment of Bridges and Wernet. 16 The streamwise velocity 
component and TKE are normalized by the fully expanded velocity and the square of the fully expanded 
velocity respectively. The spatial coordinate is normalized by the nozzle exit diameter D = 0.0508 m. The jet 
operates at Mj = 0.513 and TTR = 1.00 from the convergent SMC000 nozzle. 



Fun3D x/D=l 

Fun3D x/D =4 

Fun3D x/D=16 

PIV x/I)=l 

PIV x/D =4 

PIV x/D=16 



(a) Velocity 


(b) TKE 


Figure 6. The (a) steady RANS streamwise velocity component and (b) TKE at x/D = 1, x/D = 4, and 
x/D = 16 compared with PIV data from the experiment of Bridges and Wernet. 16 The velocity and TKE are 
normalized by the fully expanded velocity and the square of the fully expanded velocity respectively. The 
spatial coordinate is normalized by the nozzle exit diameter D = 0.0508 m. The jet operates at Mj = 0.513 and 
TTR = 1.00 from the convergent SMC000 nozzle. 
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(a) Velocity (b) TKE 


Figure 7. The (a) steady RANS streamwise velocity component and (b) TKE along the jet centerline (y/D = 0 
&; z/D = 0) compared with PIV data from the experiment of Bridges and Wernet. 16 The streamwise velocity 
component and TKE are normalized by the fully expanded velocity and the square of the fully expanded 
velocity respectively. The spatial coordinate is normalized by the nozzle exit diameter D = 0.0508 m. The jet 
operates at Mj = 1.29 and TTR = 1.00 from the convergent-divergent SMC016 nozzle. 



Fun3D x/D=l 
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(a) Velocity 


(b) TKE 


Figure 8. The (a) steady RANS streamwise velocity component and (b) TKE at x/D = 1, x/D = 4, and 
x/D = 16 compared with PIV data from the experiment of Bridges and Wernet. 16 The velocity and TKE are 
normalized by the fully expanded velocity and the square of the fully expanded velocity respectively. The 
spatial coordinate is normalized by the nozzle exit diameter D = 0.0508 m. The jet operates at Mj = 1.29 and 
TTR= 1.00 from the convergent- divergent SMC016 nozzle. 
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(a) y p /D = -10 


(b) y p /D = -8 


(c) y p /D = -6 





(d) y p /D = -4 


(e) y p /D = -2 


(f) y p /D = -1 


Figure 9. Radial variation of TKE at x/D = 10 downstream of the nozzle exit. The spatial coordinate is 
normalized by the nozzle exit diameter. Parts (a)-(f) show the steady RANS solution with the plate located 
at 10, 8, 6, 4, 2, and ID from the jet centerline respectively and extending 10 and 20D downstream from the 
nozzle ( D = 0.0508 m). The jet operates at Mj = 0.513 and TTR = 1.00 from the convergent SMC000 nozzle. 
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St 


Figure 10. The free-held prediction at R/D = 100 and ip = 90 degrees using the model of Tam and Auriault 15 
compared with the experiment of Bridges and Brown. 17 The jet operates at Mj = 1.00 and TTR = 1.00 from 
the convergent SMC000 nozzle with D = 0.0508 m. 
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(a) Xp/D = 10 


(b) y p /D = -1 


Figure 11. Attenuation plot of the free-field prediction at R/D = 100 and ip = 90 degrees using the model of 
Tam and Auriault 15 for the isolated jet and plate cases. For part (a), the plate is located at 1, 2, 4, 6, 8, and 
10 D laterally from the jet centerline and extends 10 D downstream from the nozzle. For part (b) the plate is 
located at ID laterally from the jet centerline and extends 4, 10, and 20D downstream from the nozzle. The 
jet operates at Mj = 0.513 and TTR = 1.00 from the convergent SMC000 nozzle with D = 0.0508 m. 




(a) x p /D = 10 (b) y p /D = -1 

Figure 12. Attenuation plot of the free-field prediction at R/D = 100 and ip = 90 degrees using the model of 
Tam and Auriault 15 for the isolated jet and plate cases. For part (a), the plate is located at 1, 2, 4, 6, 8, and 
10 D laterally from the jet centerline and extends 10 D downstream from the nozzle. For part (b) the plate is 
located at ID laterally from the jet centerline and extends 4, 10, and 20D downstream from the nozzle. The 
jet operates at Mj = 1.29 and TTR = 1.00 from the convergent SMC016 nozzle with D = 0.0508 m. 
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(a) x p /D = 10 (b) y p /D = -1 

Figure 13. Attenuation plot of the free-field prediction at R/D = 100 and ijj = 90 degrees using the model of 
Tam and Auriault 15 for the isolated jet and plate cases. For part (a), the plate is located at 1, 2, 4, 6, 8, and 
10 D laterally from the jet centerline and extends 10D downstream from the nozzle. For part (b) the plate is 
located at ID laterally from the jet centerline and extends 4, 10, and 20D downstream from the nozzle. The 
jet operates at Mj = 1.50 and TTR = 1.00 from the convergent SMC016 nozzle with D = 0.0508 m. 
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(a) y p /D = -10 



(b) y p /D = -6 



(c) y p /D = -4 



(d) Vv/D — —2 



Figure 14. Contours of peak acoustic source downstream from the nozzle exit at / = 1000 Hz (St « 0.3). The 
spatial coordinates are normalized by the nozzle exit diameter. Parts (a)-(e) show the acoustic source with the 
plate located at y p /D = —10, —6, —4, —2, and —1 respectively and extending 10 D downstream from the nozzle. 
The jet operates at Mj = 0.513 and TTR = 1.00 from the convergent SMC000 nozzle with D = 0.0508 m. 
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Figure 15. Results of AO AS PL as a function of non-dimensional number F (Eqn. 5). AO AS PL is calculated 
relative to the isolated jet case for each airframe configuration and jet condition as listed in Tables 1 and 2, 
respectively. 
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